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Abstract 

We consider the perturbative computation of the N-point function of chiral densities of massive free fermions 
at finite temperature within the thermofield dynamics approach. The infinite series in the mass parameter for the 
N-point functions are computed in the fermionic formulation and compared with the corresponding perturbative 
series in the interaction parameter in the bosonized thermofield formulation. Thereby we establish in thermofield 
dynamics the formal equivalence of the massive free fermion theory with the sine-Gordon thermofield model for 
a particular value of the sine-Gordon parameter. We extend the thermofield bosonization to include the massive 
Thirring model. 

1 Introduction 

The bosonization of fermions has proven in the past to be a very useful technique for solving quantum field 
theoretic models in 1+1 dimensions [1 . In a recent paper [5] we have considered the operator formulation of the 
bosonic representation of massless free fermions at finite temperature (thermofield bosonization) with the thermofield 
dynamics formalism [21 El [51 [HI HI HI [H] • The well known two-dimensional Fermion-Boson correspondences at zero 
temperature are shown to hold also at finite temperature. In Ref. [2] we have also used the thermofield bosonization 
for obtaining the real time fermion N-point functions of the massless Thirring model at finite temperature. 

The equivalence of the massive Thirring model and the sine-Gordon theory at finite temperature in the imaginary 
time formalism has been the subject of a number of authors [101 111] . The discussions have been carried out 
predominantly from the functional point of view. 

The main objective of the present work is to extend the thermofield bosonization approach presented in Ref. [2] 
to the case of massive fermions. To this end we use strictly fermionic techniques on the one hand, and bosonization 
techniques on the other, in order to demonstrate within the Thermofield Dynamics approach the equivalence of the 
theory of massive free fermions at finite temperature with the sine-Gordon thermofield theory. The demonstration 
will be done for the case of the N-point functions of chiral densities, by working in the interaction picture of the 
respective formulations, with the mass term playing the role of the interaction Hamiltonian. The selection rule 
to be imposed in both formulations (fermionic and bosonic) in order to prove the equivalence emerges here in an 
interesting way. As a byproduct the role of the "tilde" fields in thermofield dynamics and the lower branch at 
t = — i (3/2 in the real time formalism [7] emerges naturally in this calculation. 

The intimate relationship between the thermofield dynamics formalism and the algebraic formulation due to 
Haag-Hugenholtz-Winnink (HHW) of statistical mechanics has been established in Ref. [8j. In this reference the 
relevance of tilde objects to the modular conjugation appearing in the algebraic formulation of statistical mechanics 
in the HHW formalism based upon the KMS condition is clarified. This formulation of thermofield dynamics in a 
way consistent with the HHW formalism enable to extend it to gauge theories and becomes crucial in the treatment 
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of the Faddeev-Popov ghosts [5]. To this end, a "new" version for the thermofield dynamics approach for fermions 
is presented in [5]. In this paper we shall follow this revised version for the thermofield dynamics approach for 
fermions since it corrects a mistake in our previous work [5] (this is discussed in Appendix) and as we shall see, 
becomes fundamental for the success of the thermofield bosonization scheme. 

The paper is organized as follows: We begin in Section 2 by considering the thermofield dynamics approach 
to the massive free fermion theory. We use the generalization of the perturbation theory to finite temperature [9] 
in order to compute in thermofield dynamics the N-point function of chiral densities of massive free fermions as 
a power series in the mass M of the fermion, with an explicit expression for the expansion coefficients as a ratio 
of temperature dependent polynomials. In Section 3 we repeat this calculation for the bosonizcd formulation and 
extend to include the massive Thirring model. The equivalence of the two formulations is thereby established, 
upon taking suitable account of the selection rule emerging in this calculation. In Section 4 we discuss the physical 
meaning of the selection rule. We conclude in Section 5 with some comments. In Appendix we discuss some 
modifications of the thermofield dynamics formulation for fermions of Ref. [5] due to Ojima, which shall play an 
important role in order to obtain the correct bosonized expression for the Fermi thermofields. This streamlines the 
presentation of Ref. [2] . 



2 Massive Free Fermions 

In thermofield dynamics the construction of a field theory at finite temperature requires doubling the numbers of 
fields degrees of freedom by introducing "tilde" operators corresponding to each of the operators describing the 
system considered [3j [4j [SJ [6J [Tj [8l [9] . To this end, let us consider the total Lagrangian of the two-dimensional 



massive free Fermi field corresponding to the fermion doublet 



£ = £-£, (2.1) 

where £ is the usual Lagrangian of a massive free fermion, and £ is the corresponding Lagrangian in terms of the 
field iji and obtained from £ by the tilde conjugation defined by (cip) — c*?p. Since we shall consider perturbation 
theory around the massless theory, we choose to write £ in the form 

£ = £ Q + £i (2.2) 

where 

£ = ^7 M 3 M V + W'fdrf, (2.3) 

and 

£i= - - W>) ■ ( 2 - 4 ) 

At finite temperature and within the thermofield approach, the perturbative computations can be performed 
using the generalization of the Gell-Mann-Low formula for T ^ [9] ■ The vacuum expectation value of time-ordered 
products of Heisenberg operators <E> in the physical vacuum \Q{(3)) at finite temperature [71 [HI [S] is given by 



(0,/?|Tl]>fc(zfc)e X p(i / d 2 zli{x))%p) 

(oG9)|rn**(**)|o(/?)> = t , (2-5) 

k (Q,j3\Texp(i / d 2 z£ I {z))\Q,l3) 



where the right hand side is computed in the interaction picture. The thermal interaction picture vacuum is defined 

by 



1 The conventions used are: 

7 ° = (; J),7 1 = (^ 1 I) ,7 6 =7V, e 01 = l,9 00 = l,x± = x°±x\d ± =d ±d 1 . 
For massless scalar field <p(x) = <p(x + ) + rf>(x~), and for the pseudo-scalar field 4>{x) = rp(x + ) — <p(x~). 
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\0,f3) = U F [6 F ]\0,0), (2.6) 

where |0,0) is the interaction picture vacuum at zero temperature, and U f [9f] is the unitary operator (We shall 
adopt the "revised" thermofield dynamics approach for fermions as formulated in Ref. [8]) 



dp e F (b 1 i ,p) (bip 1 ) Hp 1 ) + &t (p 1 ) 6+ & ) + d ( P ) + <jt (p 1 ) <jt (p i ; 



C/ F [M = e J — V 7 , (2.7) 

where the Bogoliubov parameter 6 pip, 0) is implicitly defined by (p = Ip 1 !), 

cos 9 f (p; (3) = . 1 = , (2.8) 
V 1 + e ~ 



and 



e -/3p/2 

sin6» F (p;/3) = ___ , (2.9) 



with the Fermi-Dirac statistical weight given by, 



N F (p;(3) = sin 2 F (p;/3) = -^-y . (2.10) 



We can pass the time independent unitary operator (|2.7|) through the time ordering operation, and rewrite the 
Gell-Mann-Low formula in terms of the zero temperature vacuum |0, 0) and the transformed annihilation operators 
are given by [5] 

b(p;0) = b{p) cos e F { P ;P) + ib\p)sme F {p- 1 (3), (2.11) 

b(p;0) = b{p) cos e F { P ; (3) - itf(p)sm6 F (p;l3), (2.12) 

and their adjoints, with similar expressions for d and d. For free massless fermions in 1 + 1 dimensions the spinor 
components ipi(ip2) are left (right) moving fields: 

#"0 - ( t±-\ ) . (2-13) 



ip 2 {x ) 

and similarly for the tilde fields. Following the approach given in Ref. [S], the corrected expression for the Fermi 
thermofield is given by (see Appendix ) 



where 



^5/3) = J dp[f p {x ± ) (b( TP ) cos 9 F (p;j8) + itffrp) sin0 F (p;/?)) 

+ f;{x±) (dt(Tp) cos^b;/?) - ^d(T F ) sin 6 F (p; /3) ) } , (2.14) 

^5/3) = --= jf d F { /;(x±) (fcfrp) cos0 F (p;/3) - z6 t ( T p) sin0 F (p;/?)) 

+ (rf f (Tp) cos0 F ( F ;/3) + i d(^fp) sin0 F (p; f3) ) } , (2.15) 

= e -*f». (2.16) 
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2.1 N-point function of chiral densities from fermionic point of view 

In order to compactify the calculation it proves convenient to introduce the following notation for the Fermi fields. 
We label them by an upper index rj taking the values zero and one, with the identifications 

tJj v — if) , for i] = 0, 

iP r ' = ft ,for i] = 1. (2.17) 
Now, let us consider the chiral densities of the massive free Fermi field in terms of the spinor components 

J+i ^) = ■^\{x)^2{x): , J v jf(x) = :if>l(x)i/}i(x): , 

Jlf(x) = :Mx)^ 2 (x): , Jlf(z): = . (2.18) 

In a compact notation these expressions become 

Jl(x) = :ipfr 2 : , 3l(x) = :^f^: . (2.19) 

Here the notation for the lower index (chiral) has been chosen such as to facilitate later comparison with the 
corresponding bosonized expressions. In terms of the chiral densities the mass term reads, 

m$h> - w>) = m E j a • ( 2 - 2 °) 

A=±l »j=0,l 

Note that the minus sign on the left hand side for the tilded fields has been taken into account by reordering the 
fields on the rhs, taking account of the definitions (|2.17j) and of Fermi statistics. With the identification l|2.4[) for 
the interaction Lagrangian we have, setting J = J , 

(0(/3)|TJ +1 (x 1 )...J +1 (x A r)J_ 1 (y 1 )...J_ 1 (y w ,)|0(/3)) = 
j^T0) E n! E (O,O\TJ+i(x 1 ,p)...J +1 (x N ,{3)J- 1 ty 1 ,0)...J- 1 (y N ,,p) I dZ n J[ ( J^ k (z k ,^)) \0,0) , 

F "=0 ' (Afc,%) ' fc 

(2.21) 

where we transfcrcd the temperature dependence of the ground state to the densities by commuting the time 
independent unitary operator Uf[0f] through the time ordering operation. J±i(x; f3) are the thermal chiral densities 
of the free massless Fermi field 

J+i(x;0) = UF^J+iix^iep] = ■4\{x + -f3)4, 2 (x--f3):, (2.22) 

J-i(x;p) = UF[eF]J-i{x)Up X [e F ] = itl>t(x-;p)ih.(x + ;p):, (2.23) 
the integral J dZ n is short hand for 

/ n />oo 
dz n = n / d2zk ' ( 2 - 24 ) 
fc=1 J -co 

and N"f(P) represents the contribution of the vacuum graphs, 

WHE^r"- E /^ n (0,6|rn(j^(^,/3))|0,6). (2.25) 

We can thus apply Wick's theorem to compute each term in the expansion in powers of M. The two-point functions 
for the fermion thermofield have been computed in Ref.[5] and adopting the notation given in (|2.17[) the fermionic 
propagators are given by 
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<0,0|T^WV;/3)|0,0) = 2lp ^ [¥x± _l ± _ l£<x0 _ y » ))V (2-26) 



and 



(0,0|T^°(a : ± ;)8)V lt (l/ ± ;)9)|0,0) = - (0, O^ 1 ^; W V; /3)|0, 0) = — cosh[ ^ ± _ y±)] ■ (2.28) 

These functions can be collected in a compact notation by using 

cosh(x) = ±i sinh(x T «7r/ 2 ) (2-29) 

so that we can write 

(o, o|iy (xf ; m m 1 (yf ;0) |o, o) = i— r- — — -tt • ( 2 - 30 ) 



2i (3 sinh 



/3 ri 



In applying the Wick's theorem only terms with equal number of ip^ and ip^ survive irrespective of their -q 
values. The same is true for the second spinor component. In terms of the chiral densities this leads to the 
selection rule requiring the sum of all values of A have to vanish, irrespective of the r\ values. There are thus equal 
number of positive and negative values of A. We shall denote the space-time coordinates of the fields associated to 
A = 1 by Xi , while the ones associated to A = — 1 are denoted by yj . The values of the r\ upperscript are accordingly 
splited as rji and rjj . The result of the computation can be written as 

n n 

(0,6|T JJ Jl =+1 ( Xi ; 0) JJ 4 , • //, ; (3) |0, 0) = 

i=i j=x 

^ - irnP/2) - (Vj - *hP/W) - irnP/2) - (y+ - ' (2 ' 31) 

where 

lo( Xi - Vj ) = (2*/3) sinh-^ - x.j - i(~l)^ee(x° t - y°)) . (2.32) 

Note that the factors (i) ni+rij can be factorized out of the determinant computation resulting in an overall sign 
depending on the total numbers of tilde fields if>i and of corresponding tilde fields ip2- Both contributions lead to a 
factor (-l) m , with 



m 



representing the total number of J± in the string. 

Now using the factorization formula [12] for the determinant of the n x n matrix l/tv(xi — yj) 

JJ uj(x 1 -xe) JJ u(yj - y r ) 
det = ^ j<j ' , (2.34) 

and suppressing for a while the it prescription, we obtain 
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(o,oirn^ 1 (^)n j ^^)i°'°)= 

i=l j=l 

n n 

JJ cj((4 - %/3/2) - (x+ - irH'0/2)) U ~ Wrf/I) - (Vj> ~ %0/2)) 

j<j' 

71 



j<j' _ 

n 

n n 

n - - - ^OT) n - w 2 ) - - ^-'Z 3 / 2 )) 

j<j v 

n 

n^((^-^/3/2)-(%-*^/3/2)) 



(-l) m — s — . (2.35) 



Here we introduced the function 

= wfc+^fc - ) , (2.36) 
Returning the ie prescription means that in the denominator of eq. (|2.35p the change should be performed 

sinh(^(a;+ - yf)) 8inh(~(x7 - yj)) — > sinh(^(x+ - yf - ie^{x\ - y°))) sinh(^(a;7 - yj ~ ie^x? - yf))) , (2.37) 
where Ej = (—l) nj e. 

Returning now to the expression (|2.21|) we reorganize the terms in the expansion (|2.21|) according to chirality, 
with n being the total number of internal currents with A/- = 1 and n' the same for A; = — 1. To the former we 
associate the variables Zk and to the last ones z[. The contributions can now be collected into the expression 



N N' 

>n 



(o(/3)ir(nj + i(^)n j -ife-))i°^)) 

i 



n+n'=0 (Vk,vi) 



(O,O\TJ +1 (x 1 ,0)..J +1 ix N ,0)J-x(vi,P)...J-i{vw (2.38) 

k I 

where the notation * means that the sum obeys the chiral conservation condition that now reads 

N — N' + n — n' = 0, (2.39) 

and dZ n refers to integration over the variables {z 1; • • • , z n }. Each term of this expansion is obtained directly from 
(gHHD leading to 



N N' 

(0(/3)|r(nJ + i(^)Il J -ife))|0(/3)) 

i=l 3 = 1 
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E* ( E / d W dZn ' (-!)*!(*, K*,/), (2.40) 



^ n+n'=0 ( % ,^) 



with 



I(x,y;z,z') 



N JV' 

i<i' = l j<j'=l 
JV JV' 

UllSlixi-yj) 

i=ij=i 



JV n JV' n' 

n n - - *^/3/2)) n - ( z '< - *^/ 2 )) 

t=l fc=t 3 = 1 ' = 1 

TV n n JV' 

nn^ - [z\ - iprn/2)) n n ^ - ( Zk - m?^/^ 

i=li=l fc=l j=l 

II n((zk-iVkP/2)-(zk>-irik>0/2)) J[ Q((z'i - i m (3/2) - (z' v - i m ,{3/2)) 

k<k> = l Kl' = l 



(2.41) 



AT JV' 

fc=i i=i 

It is interesting to note, that the shift in the argument of the two-point function involving tilde field, (|2.30p , 
(|2.3ip . can be understood in the context of the real time formalism as the tilde field living on the lower branch 
of the integration contour localized at Im(t) — — z/3/2 in the complex time plane. This feature appears further in 
expression (|2.41[) since to all integrated variables z k are associated fields with r] labels and these labels are summed 
for r\ — and rj = 1. In the real time formalism this corresponds to contributions from the upper branch, the lower 
branch and functions from mixed branches. This has been observed to be necessary to satisfy the KMS conditions 
®- 

Note that all expressions in the numerator link currents of the same chirality, whereas all terms in the denom- 
inator link currents of opposite chirality. We thus associate with each coordinate £j the chirality A^. We further 
reorganize the terms within the products by separating the integration variables into those which belong to ordi- 
nary currents and tilde currents (i.e. z k and z k , respectively), and denote the respective chiralities by Xi and — \ k . 
The assignment of — \ k to the tilde currents chirality is due to the definition Eq. (|2.18p . With this we can write 
expression (|2.4ip in the following form, suitable for later comparison: 

JV JV' oo 

(0G8)iT(nj+ifo) n j -ife)) = TT7^E(- iM ) nx 

i=l j=l n=0 



E^if(-ir/n E*^/IP 

m,m 1=1 {A,,A fc } fc=l 



' 2 Zk x 



JV JV' 

1] [nfc-a*;/?)] J] [n( yj - yjl ;f3)] 

i'>i f>3 
N N' 

n n [n(xi- yj ;i3)] 
%=i j=i 

JV m JV m .q 

n n - *« Ai n n - ** - f )] _ Afc >< 

i=l £=1 i=lfe=l 
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JV' m AT' m 'fl - 

n n - p)\ A n n Hvi * f )] Xk >< 

i=i £=i j=i fc=i 

where 

fifa - x,;/?) = (2z/3) 2 sinh [~(x+ - as+ - i £ e(x? - a:°))] sinh [-(sT - X J ~ ise(x\ - x^)\ , (2.43) 



(x, - Xj ;0) = (2z/3) 2 sinh [^(x+ - xf + i£ e(x° - x°))] sinh [^(x7 - x7 + J£e (x° - x°))] , (2.44) 

Sl(xi - Xi - f ) = (^) 2 sinh [^(x+ - 4 - |)] sinh [Ifa xJ - |)] , (2.45) 
and the sum over the chiralitics Xg and respects the total chirality condition 

N - N' + J2 x e-J2~ Xk = - ( 2 - 46 ) 

£ k 

This concludes the computation of the N-point function of the chiral densities in the fermionic version. In the 
next Section we shall consider the computation of these N-point function within the bosonized version of the theory. 

3 Thermofield Bosonization point of view 

The next step is to compute the N-point function of the chiral densities of the massive free fermion theory from the 
thermofield bosonization point of view. Since at T = the bosonized theory describing both the massive free fermion 
theory and the massive Thirring model is the sine-Gordon theory |T3l [14] with distinct values of the sine-Gordon 
parameter k ( different values of the scale dimension of the mass operator ), we shall consider the perturbative the 
perturbative computation of the N-point function of chiral densities of the bosonized massive Thirring model from 
thermofield dynamics point of view. 

To begin with, let us consider the sine-Gordon theory for the doublet (<£>, $), whose Lagrangian can be decom- 
posed as 



with 



£(x) = £ (0) (x) + £ 7 (x), (3.1) 

£ (0) (x) = C w {x) - C w (x) = ^:d^Mx)d^<S>{x): - 1 :d^(x) d»®(x) : , (3.2) 

Zj{x) = Af (-) (:cos(/c$(x)): - : cos (k $(x)) :) , (3.3) 

where fj, is the infrared regulator (IR) reminiscent of the free massless scalar theory. The scale dimension D of the 
mass operator is given in terms of the Thirring coupling parameter g as [13] 

D-±. (3.4) 

For the values k 2 — An, the sine-Gordon theory describes the massive free fermion theory discussed in the preceding 
section. 

Now, let us introduce the Mandelstam representation [TBI ITS] for the Fermi field operator 
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POO 

*(z) = e 4 ' :e Jx* :, (3.5) 



/•CO 

= e 4 :e ^i 1 : . 



(3.6) 



Notice that the field operator VP is not obtained from ^ by the "tilde" conjugation operation. This follows from 
the fact that ^ should be an identical copy of <J/ carrying the same charge and chirality quantum numbers. This is 
explained in detail in Section 4 and in the Appendix. 
The bosonized chiral densities J±i(x) are given by 

J+i(*) = -©WW, (3.7) 
J-iW = -(£)w*(z), (3.8) 



,2vr 

where the Wick-ordered exponentials carrying opposite chirality are 

W(x) = :e 4K * (x) :. (3.9) 

W*(x) = :e- iK *^:. (3.10) 
The interaction picture vacuum is now given by 

\o,P) = u B [e B ]\o,o), (3.ii) 

where the unitary operator taking one to the thermofields is given by dp 1 ] = p) 

U B [0 B ] = e" JT ^(^>^>-*V)«V))#.W) f (3 12) 
and the Bogoliubov parameter 6 B (p,/3) is implicitly defined by 

_ 0p 

sinh0 B (p,/?) = 62 (3.13) 
Vl-e-^P 

cosh0 B (p,/3) = - 7 =L==, (3.14) 
with the Bose-Einstein statistical weight given by 

N B (p,(3) = sinh 2 B (p,/3) = ^—j- (3.15) 

Following the same procedure of the preceding section, in the bosonized theory the N-point function of the chiral 
densities is given by 

N N' 



<oos)|r(lJj+i(^i) II J -i^)) l°C9)> = 
»=i j=i 

N N' 

(0,0\T(j[j +1 ( Xi ;(3) H J_i(y j; /3) exp / £,(z; /?) d 2 z}) |0, 0) 



i=i i=i 



(0,0|Tcxp{i y £/(z;/3)d 2 z||0,0) 
where the interaction Lagrangian (in the interaction picture ) at finite temperature is 



(3.16) 
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C t (z;P) = UbWi^CiWUb^Ob] = M (J- e - 2Z? ^')) (: cos ( K 0(z; /?)) : - : cos (k^z; /?)):) , (3.17) 
J±i(x; (3) are the thermal chiral densities in the interaction picture [2] 

J+i(xi;0) = U B [eB]J + i{xi)Uz X [9 B ] = - (^e~ 2D z ^)w{ Xi] ^\ Xi ) , (3.18) 

J_!( yj -;/3) = [/bNL!^)^ 1 ^] = - (^e- 2 ^^')) H" 8, A,, i. (3.19) 
with the thermal Wick-ordered exponentials carrying opposite chirality given by ( in our convention X Xi — X yj = 1) 

W(xi;/3,\ Xi ) = :e < ' eA »* (3.20) 

W*(y 3 ;f3,\ y] ) = :e - 4 « A «^fc-/3) :j (3.21) 
<j)(x; (3) is a free massless pseudo-scalar thcrmofield [2] and z(/3, //) is the infrared divergent integral [2] 

Expanding the exponential in powers of M, introducing the tilde Wick-ordered exponential 

W?*(z;/?,A 2 ) = (3.23) 
and using the fact that the fields <j)(x;(3) and 4>{y,0) commute [2], we get 

i J d 2 z C(z) i M (j- e~ 2 D z ) { J d 2 z: cos k0(z; /?) : - J d 2 z : cos «0(z; /?) : | 
(i^(^ e - 2Z5z )"( / d 2 z: cos^(z;/3): - f d 2 ~z : cos K 0(z; 0) : 

n=0 ^ *' 77 J J 



{tM)n (^e- 2 ^)" g ( " ! ^ B (~ir ( / cos^(z;/3)=) m ( J d 2 ~z = cos K fe : )' 



E (n!) 



m.m 



oo „m„mm m 

n=0 " m ~ ' ' ^ f=l ^ fe=l 1=1 k=l 



oo j. „ m „ m 

n=0 ~ J *=1 17 fc=l 



2 Zfc X 



{M™ {A*}™ 1=1 fc=l 

where Af , A& = ±1, and ( ) runs over all possibilities in the set {Ai, . . . , A m }({Ai, . . . , A~}). Denoting 

{Mm {A fc }„~, 

by Nb{(3) the normalization factor, the Green's function (13.16|) can be written as 
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N N' 

<o(/3)ir(nj+i(si) n j -i(%)) w)) 

i=l j=l 



nhi tt-r^MT (ii e — E^ff (-!)"■ / ft / ft * 



2 z k x 



n=0 ~ 1=1 k=l 



N N' mm 

J2 (O^T^W^;/?^ ]JW*(y,0;^^ (3.25) 
The time-ordered product of two Wick-ordered exponentials is denned by, 

: Wfo /?, X x ) W(y; 0, X y ) : ( (w(x; f3, X x )W{y; 0, X y ))o(x° - y°) + (w(y; /?, X y )W(x; (3, X x ))9(y° - x°f 

= : W(ar; /3, A,) W(y; (3, X y ) : (t(W(x; (3, X x )W(y; /?, X y )) ) , (3.26) 

where 

(w(x;(3,X x )W(y;(3,X y )) = e~ x * x * K * (o,o|*(*;/3)*(»i/J)(o,o> _ (3 27) 

Using the identity (0(ac 4 ) = 6>(a;° - y°) = 0y) 

e (<f>i<f>i)() i . _|_ e^^Oji = e^^^ 9 ^ + ^ j ^^ 9ji = e^^^^ (3 28) 

we can write 



T 



(w(x; [3, X x ) W(y; (3, X y )) = : W(x; (3, X x ) W(y; 0, X y ) : e~ x ^ ^ (T*toMv,f») . ( 3 .29) 

The Wick's theorem can be extended to generalized time-ordered product of Wick-ordered exponentials and we 
obtain 



T(J[W(xj;/3,Xj)]lW(x k ;l3,Xk) 

3=1 k=l 



[] (T{W{xf,t3,Xj)W{x r -,i3,X r ))) \{ (T{W{x k -,(3,X k )W{x k r, l3,X k '))) Y[(T(W(x fl (3, X 3 )W(x k ; (3,X k )) f 

j>j' k>k' j,k 

(3.30) 

In this way one gets for the Fock vacuum expectation value of T-ordered product of Wick exponentials in (|3.25p , 

N N' 



{Q,O\T(j[W(xi;0,Xi) n w *(%;AAi) [] W(z e ;(3,Xt) J[ W*{z k \ P, A*)) |0,0) 

i=l j=l 1=1 k=l 

N N' 

i I (T{W{x i ;p,X i )W{x il ;P,X i ,))\ J[ (t (W* (yj ; (3, Xj)W* (yj> ; ft, Xj> )) * 



i'>i 3>J 
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N N' 

II II (r(^(^;/3,A 4 )^*( % ;/3,A,)) 
i=i j=i 



N m N m 



n n (T(w(xi;p, Xi)w(z r , p, ) n n & p> **)) 



2=1 1=\ i=l k=l 



N m N m 

m rn 

Y[ (T(W(z i ;l3,Xe)W{z il ;p,X if ))) J| (r(W*(zk; P, Xk)W*(zk r ,P, Xk')) 

£'>£ k'>k 
rn m 

nil( r (^^>#feAA*)))- (3-31) 

The propagators of the scalar thermofields can be written in terms of the propagators of the Fermi thermofields as 
follows g] 

(O,Q\T4>(x;0)4>{ V ]P)\O,O) = /3) - ±- In (£) - -L lnJl(x - y;/3) , (3.32) 

7T Z7T Z7T 47T 

(O,O|T0(x;/3)^;/?)|O,O) = -i + ^z(ji',/3) - ^M^) - ^ lnfi(z-y;/3) (3.33) 

(O,O|T0(x;/3)^(y;/3)|O,O> = - -L /(//', /3) - i- In (2/3) + i- lnft^ - y - ^) , (3.34) 
where the f2's arc defined in (|2.43p - (|2.45p . The dependence on the infrared cut-offs // and fi" are given by [5] 

z(»',P) = [°° *(e* - l)" 1 , (3.35) 
P 



oo 



dp 1 
7^rnh(f ) 



/V'>/3) = / --T7ST' (3 - 36) 



and the corresponding asymptotic behavior are 



*(j*'«0,j8)- -L + ihW), (3.37) 
PM 2 

/V'wO.fl--^. (3.38) 
The N-point function is then given by (A ai = A % = 1) 

N N' oo 

(0C9)|T(nJ+i(*i) Il J -i^)) |0(/3)) = ^vgr £ (-l) iN+N,) (iM) n x 



<=i i=i 



E^f(- i r E E F^M,»»,p)G^M,ff) / / n^ x 



=1 " fe=l 
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where the phase 



N 



N' 



n m & 



%'>% 



j'>3 



N N' 



n n Mxi-wp)] 



DXiXj 



= 1 3 = 1 



N 



n n - ^^)] DXiXi n n ^ - - t>] ~" AsAfc x 



AT 



i=l «=1 



JV' m 



DXjXi 



i=l fe=l 
AT' m 



nn 

i=i fc=i 



nn 

m m 

[] [fi(z,-^;/?)] DA ' V [] [Q(~z k -~z k >;P)] 



«/3.-| DXjX k 



DX k X k , 



e>>e 



k'>k 



K. 



mm ■ n ~ 

n n [n^ - ^ - f )] m ^ 



i k 



2iirD AfeAfe' 
X = e fe '> fe 

is the identity for integer values of the scale dimension D and the cut-off dependence is given by 

m m 

d(n -n' + E a ^-E a *) 



(1 - D)(N+N'+r, 



2tt 



e=i k =i 



m m 

(VA fe )(iV - N' + VA,) 



) 1 

v 7T fl' 



N - N' + 



In arriving at these expressions we collected the cut-off dependent terms, noting that 

N N' N N' N m JV' m m m 

e + e ^j^j' - e X{ e x 3 + e a * e a ^ ~ e a j e Xe + e a ^ a£ ' + e AfeAfe ' 

i'>i j'>j i 3 it, j I £'>£ k'>k 

-(jv - TV + E A ^-E Afc ) - 2^ +7V ' +n ) + Y,~ Xk { N - N> + E\) > 



=i *;=i 



with (Aj = Aj : = 1) 



JV 



E A * = ^ = E( A *) 2 ' 



i=l 
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N' N' 

5>= N' = J2(*j)\ (3-45) 

and made use in (|3~f2f of the asymptotic behavior (|3~37f and (pT38j) . 

Now, let us consider the free massive fermion theory (D — 1), which is an infrared cut-off independent scale 
non-invariant theory 

Following the procedure introduced in Ref. [5] (see Appendix), the theory of the free massless scalar thermoficlds 
<p(x; P) and 4>(x; j3) can be constructed as the zero mass limit of the massive free scalar thermofields E(cc;/3) and 
S(x; [3). In this way, the infrared regulator fi of the zero temperature two-point function should be identified with 
the infrared cut-offs \i! and /i" of the temperature-dependent contributions z(//,/3) and /(//', /3), 

M " = = - , (3.46) 

such that 

i^™WV',/3) =1. (3.47) 

M' — A* — £ 

In this way, the only non zero contributions in the expansion (|3.39p are those which satisfy the selection rule 

?71 771 

N - N ' + H Xe ~ ~ Xk = ■ ( 3 - 48 ) 
t=i k=i 

The selection rule (|3.48|) is in agreement with our previous computations (Eq. (|2.39jl ) for the massive free 
fermionic theory (D = 1). Indeed, for D = 1 the bosonized version of the perturbative expansion given by (|3.39[) 
coincide with our previous result in the fermionic formulation obtained in Section 2. In order to recover the factor 
(— 1)" that appears in (|2.42[) . we shall make use of the selection rule (|3.48p . Taking into account the summations 
in (|3.39[) . one can write 

{Mm { Afc 

(-1) ( iM )« = (-l)( Ar+A "+")(- l M)" ee (-^"^(-iM)" ee 



mm mm 



(N - N' + x e + Xk ) (N-N' + J2 X e ~Yl Xk ) 

(-1) t=i fc=i (-iM) n =(-l) *=i fc=i (-iM) n = (-iM) n . (3.49) 



In this way, for D = 1, the bosonized iV-point function (|3.39p can be written as 

N N' 



<o(/3)|r(nJ + i(*i) Il J -i(%)) iocs)) 



i=l 3=1 



a* Ew'"^ £ £ /n^/lK* 



Zfe X 

=1 " fc=l 



AT JV' 

IJ [n{xi-Xi>;0)] TJ [fi( % - % ,;/3)] 
j'>j 

N N' 

n n [n(xi-y i; /3)] 

»=i i=i 
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N m N m 

n n ze; /?)] Xi n n - ^ - 1 r Afc * 

i=l i=l i=l fc=l 

AT 7 m TV' rn . n ~ 

n n t n ^- - ^ A< n n ^ -^-t^ Xkx 

j=ie=i j=ik=i 

TfL _ _ (It fit . n _ 

\\&\ei i r rPw~ ~ /oxi^fc-^fc' TT TT fr^/ - — 



n - z e; 0)]^ n [n(i fe - % ;/ 3)] Afcv n n - - f)] - ( 3 - s °) 

in agreement with the corresponding function (|2.42p obtained in the original fermionic version of the theory. This 
proves the formal equivalence of the massive free fermion thermofield theory with the sine-Gordon thermofield model 
for the particular value of the sine-Gordon parameter k 2 = Air. 

4 Physical Meaning of the Selection Rule 

In this section we shall discuss the physical interpretation of the bosonized expression (|3.6|) for the field *S? and 
its connection with the interpretation of the selection rule p.48j) . To begin with, let us consider the massless 

free fermionic theory described by the doublet ( .y* ] . The corresponding bosonized expressions are given by 



p.5p ~ p.6p with k = 2-^7?, and can be written as 

= (£) i :e ^{*)+^)} ;] (4.1) 

ifo) = (Jty :e ^{7 5 ^) + $(x)}. t (4.2) 
with e^d^ip = dfj,(f>. The bosonic doublet (</>, <fi) is described by the total Lagrangian 

£= ^(M* \( d ^f ■ (4-3) 

The physical interpretation for the fact that the bosonized expression (|4.2p for the field ip is not obtained from the 
corresponding bosonized expression (|4.ip for the field ip by the "tilde conjugation operation" is the following: In 
the Thermofield Dynamics formalism, the fictitious "tilde" system should be an identical copy of the system under 
consideration, which implies that the field ip should be an identical copy of ip, i. e., carrying the same charge and 
chirality quantum numbers. At T = 0, we obtain the following equal-time commutation relations for the free scalar 
fields! 

2 The corresponding canonical momenta are 

n(x) = d <f>(x) , 

U(x) = - d <t>(x) , 

in such a way that the canonical equal-time commutation relations are given by 

[<f>(x) , r%)] = iSix 1 - y 1 ), 

[?(*), 5(y)] = iSix 1 -y 1 ). 

The dynamical equations are 

da<t> = - i [<l>(x) ,H] = - i [<l>(x) , H] , 
do4> = i \4>{x) ,5] = - i [00) , h] , 

where the total Hamiltonian H = H — H is the generator of time evolution of the combined system. 
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[4>(x) , dafty)] = i Six 1 - y 1 ), 



(4.4) 



$(x),do$(v)] = -i Six 1 - y 1 ). (4.5) 
Within the Thermofield Dynamics approach the above algebraic relations are retained at finite temperature. The 

computation of the bosonized expression for the fermionic currents = tpj^tp and = ip^ip of the free massless 
theory were performed in Ref. [3J. Using the bosonized expressions for the free massless fermion fields (|4.1[1 - 
one finds 





Jfi(X,P) — — 


-j= d^ipix] 0) , 

A 


(A R\ 




Jn (X- 0) = + - 


■ — 3,.u5( x' 

i — ^ u,^ y-^i h* } ■ 


(4 7) 


The axial-vector currents are (^ 


— 6 lv) 








J 5 (x:0) = - - 


— d,,d>( x' 0) 


(4.8) 




J^\ x i0) — + - 


—j= Ofj.<P(X,p) . 


(4.9) 


Introducing the corresponding char^ 


;es 






Q(0) = 


/■+00 

/ Jniz;0)dz 1 , 

J —oo 


Q(0) = / Jn(z:0)dz 1 , 

J —oo 


(4.10) 


Q 5 (0) = 


/•+00 

/ 4iz-0)dz 1 , 

J — oo 


/"+00 

Q 5 (/3) = / Jl(z;0)dz 1 , 


(4.11) 


and using l|4.4p-(|4.5p we get 










[Q(0),<l>(x;0)] 


= -i>{x;0) , 


(4.12) 




[Q 5 i0),iPix;0)] 


= -7 5 V>(z;/3), 


(4.13) 




[Qi0)Jix;0)] 


= - ^(z; /3) , 


(4.14) 




[Q 5 i0)Jix;0)] 


= -f^x-0), 


(4.15) 



implying that ij) and ip carry the same charge and chirality quantum numbers. The total charge operators corre- 
sponding to the fermionic doublet of the combined system are given by 



Q(0) = Q(0) - Q{0) , 

Q 5 (0) = Q 5 (0) - Q 5 (0). 
For the fermionic charge of the combined system we obtain the following selection rules 

n h n h 

[Q(0), T[^0)T[^0)] =~{n~h) Y[^(x i; 0)l[if(y 3 ;0) 7 

n h n h 

[Q(0), Y[i>(xf,0)l[ft(yf,0)] = -(» + n) H^(x l ;0)Y[^(y J ;0). 

i—1 _7=1 l — l J — l 



4.16) 
4.17) 

4.18) 
4.19) 
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This implies that in the massless free fermionic theory one finds the following off-diagonal 2n-point functions 



(o,oinv'K ± ;/3)n ?/, ^ ± ;' 3 )i o ' ) 



n n 

[n(2*/3)smh^(4 - xf)] [n(2^)sinh|(y± - yf) 



f<3 



i=l 



R 

[U(2i0) S mh^(xt -yf -if -ie) 



and 



(o,oinv'(^;/3)n^ t ^^)i o ' ) = °- 

i=i j=i 

For the axial charge we obtain the selection rules (here J = ip\ip2 — :e 21 ^^: , J — xl>\ip2 = :e 2l ^^i) 

n h n h 

[S 5 (/3), l[j(^,P)l[Ayj,P)] = -2(n - fi)Y[j( Xl ,P)l[j(y 37 f3), 



(4.20) 



(4.21) 



(4.22) 
(4.23) 



Z— 1 7 — 1 4 — 1 J — 1 

This provides a clear understanding of the physical meaning of the selection rule (|3.48[) , that is, the only non zero 
contributions in the expansion (|3.39[) are those with zero total chirality Q 5 (/3) 



m m 



fe=l 
Q 5 



(4.24) 



of the combined system. 



5 Concluding Remarks 

The main objective of this paper was to prove in thermofield dynamics the equivalence of the theory of massive free 
Fermi fields to the the sine-Gordon theory for a particular value of the sine-Gordon parameter k 2 — An. Approaches 
of other authors [1Q[ 112] differ from ours in two respects: i) they make use of the imaginary time formalism, and 
ii) treat the fermionic side in a hybrid way. We have treated the fermionic side of the problem strictly from the 
fermionic point of view. On the bosonic side we have used the thermofield bosonization [2] in order to compute the 
n-point functions of the thermal chiral densities. On the bosonic side we first obtained the n-point functions of the 
chiral densities from a generalized Mandelstam operator ( with non-canonical scale dimension ) for the corresponding 
bosonized expressions at finite temperature, and then recovered from there the corresponding n-point functions of 
the free theory as a limiting case. A gratifying byproduct of our analysis was the observation, that these n-point 
functions showed in a natural way, that the tilde fields of the thermofield dynamics could be regarded as living on 
the lower branch of the integration contour in the complex time plane, displaced from the real time axis by — iw, 
in accordance with the work of ref. [7] . 

We recognize that the "revised version" of thermofield dynamics formulation for fermions introduced in Ref. 
[8] enables to obtain the correct bosonized expression for the field ip and becomes crucial in order to establish the 
two-dimensional Fermion-Boson mapping within the thermofield dynamics approach, as well as, to obtain a clear 
understanding of the chiral selection rule of the combined system. 

Acknowledgments: We are grateful to Brazilian Research Council (CNPq) for partial financial support and 
to the FAPERJ(E-26/170.949/2005)-DAAD scientific exchange program which make this collaboration possible. 
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Appendix 

Thermofield Bosonization of the Free Massless Fermion Field Revised 

In this Appendix we shall consider the revised thermofield dynamics approach for fermions presented in Ref. [8] in 
order to correct a mistake in our previous paper [2] and give the revised thermofield bosonization prescription for the 
free massless Fermi field doublet (ip, —iijj'). Although in the fermionic formulation the use of the new version does not 
change the computation of the diagonal two-point functions of the free massless fields (0, 0\ip(x ± ; 0)ip^{y ± ; /3)|0, 0) 
and (0, 0\ip(x ± ; (3)^ {y ± \ /3)|0, 0), it corrects a "sign" in our previous computation [2] of the off-diagonal contribution 
(0, 0|i?A(a: ± ; f3)ip(y ± ', P)\0, 0), and as we shall see, gives a new insight into the thermofield bosonization scheme. 

In the "old version" , the operator taking one to the Fermion thermofields is given by 03 [H 03 [5] 



U f (9f) = exp{ - / d P F (\p 1 \,f3)(b(p 1 )b(p 1 ) - &V)&V) + dip 1 ) d(p) - S (p 1 ) <fi (p 1 ))} , (A.l) 

and the corresponding transformed annihilation operators are given by 

b(p;[3) = b(p) cos 6 F (p; (3) - P (p) sin 9 F (p; 0) , (A.2) 

bip; 13) = b(p) cos F (p; (3) + & f (p) sin 6 F (p; 0) , ( A.3) 

with similar expressions for d and d. As stressed in Ref. [8], the requirement for (|A.2j) and (|A.3j) to be consistent 
with each other implies the following "tilde substitution rule" for fermions (b = b(p)) 

b = - b. (A.4) 

In the "new version" [8], the vacuum state |0(/3)) is obtained from the Fock vacuum |0, 0) by the "modified" unitary 
operator U F (6 F ), which can be formally obtained from the old one (jA.lj) by the substitution 

b^ib, (A.5) 

(A.6) 

with similar substitution for d, i. e., 

U F {6 F ) = exp{-i y°° dpepdp'lftfap^bip^+b^p^ip 1 ) + d{p 1 )d{p)+d\p 1 )d\p 1 ))). (A.7) 

The revised fermionic transformations are given by [8] 

6(p; (3) = bip) cos8 F (p; (3) + iV(p) sm8 F (p; 0) , (A.8) 

b(p;f3) = b{p) cos 9 F (p;f3) - i b\p) sin F (p; (3) , (A.9) 

with similar expressions for d and d This procedure replace the "tilde substitution rule" (|A.4[) by the "tilde 
conjugation operation" 

Jib) = -ib, (A.10) 

analogous to the case of boson operators. Taking this into account, the revised expression for the two-dimensional 
free massless Fermi thermofields are given by 

^■(3) = -L f dpif^) (b(Tp) cosO F {p-!3) + iV(Tp) sin0 F (p;(3)) + 

(cP(tp) cos F (p;0) ~ idiTp) sin^(p;/3)) } , (A. 11) 
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<Kx±;/3) = -L r dp fop) cosMP;/?) ~ ^(Tp) sin0 F (p; /?) ) + 

/p^) cos f (p;/3) + id( TP ) sin0 F (p;/3)) } , 



(A.12) 



where 



/ p (z) = e-^. (A.13) 

Now, let us compute the off-diagonal two-point function within the fermionic version. Using (|A.11[) and (|A.12[) , 
one gets[f| 



1 

2^ Jo 



dp 



(o,o|iV(^ ± ;/3)V'(y ± ;/3)|o,o) = 

which can be written as 

(O,O|i^(a: ± ;/3)V(tr fc ;0)ftO) 



cosp(:r — y^) 



1 



cosh 



2/3 cosh 5(x± - y ± ) ' 



2i/3sinhjg(x ± - y± - if) 



(A. 14) 



(A.15) 



In order to establish the Fermion-Boson mapping, let us compute the two-point function above from the bosonized 
point of view. To this end, we define 



\ 2n / 



(A.16) 



(A.17) 



where Kp is a Klein factor |16j which ensures normal anticommutativity between tilde fermion thermofield com- 
ponents and non-tilde ones [2], and the value of 7 = ±1 will be fixed at the end of the calculation. Using that 

m 



(0,01^5/3)0(^^)10,0) = -i- /(//',$ + i-ln{ cosh^a* - y*)]} , 
we obtain (the global minus sign arises from the Klein factor) 



(0,0|z^(x ± ;/3)^(y ± ;/?)|0,0) = - — e 



= _ ii„-2j ! ( M , ,/3)+ 7 /(/ 1 ",/3) 



cosh-(a; ± - y ± ) 



Using the asymptotic behavior 



*(A/»0,/3)-> + i ln(/V): 
fit*" -0,0) 



we get 



(0,0|^(;r;W(y;/3)|0,0) 



v7r/i' 



2/3 [ cosh 5 (x± - y 



±M 7 



(A.18) 

(A.19) 

(A.20) 
(A.21) 

(A.22) 



3 The Eq. (4.23) of Ref. [2] is correct, but it cannot be obtained from (4.14) and (4.13) (where we have used the old version with 
| |A.2| | and l)A.3| l'). which gives 



sinp(x ± — y ± )Np(/3,p) e 2 dp 



lip 



instead of cosp(x ± — y ± ) in the integrand. 
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In order to recover the same space-time dependence as in (|A.14[) we must require that 



7 = 1. (A.23) 

In order to obtain an infrared cut-off independent two-point function, the free massless scalar thermofield theory 
should be considered as the zero mass limit of the massive free scalar thermofield 

(0,0|E(s; J 8)E(y;/3)|0,0> BW o - (0,0|^;/3)^/;/?)|0,0) , (A.24) 
(O,O|S(s;/3)E(tf;0)|O,O> BwO - <0, O|0(ar; y9)0(y; y9)|0, 0) , (A.25) 

(0,0|£(i;/3)E(y;/3)|0 ) 0> wl ^ — (0, O|0(x; /3) |0, 0) . (A.26) 
In this way, the infrared regulator fj, of the zero temperature two-point function should be identified with the infrared 
cut-offs jj! and fi" of the temperature-dependent contributions and /(//',/?), i. e., 



m" = m' = -, (A.27) 

7T 



and we get, 



in accordance with (|A.14[) . Since the off-diagonal selection rule carried by the Wick-ordered exponential requires 
that 7 = 1, one concludes that the bosonized expression for ip(x) is not obtained just by the tilde conjugation 
operation (W) of the corresponding Wick-ordered exponential W defining ip{x). Besides a multiplicative factor 
(— i) and Klein factors, the bosonized version of the field ip(x) is obtained only by the tilde conjugation of the 
creation and annihilation components of the field (j>[x) at the exponent, which can be achieved by defining ip in 
terms of the Wick exponential W* , 



i#r;/3) = U B (6 B )iiP(x)U£ 1 (6 B ), (A.29) 



with 



,D 



i$( x ) = (fL_Y KW *( X ) = (!L-Y K . e i2^4>{x)._ (A30) 

V 2tt 1 \ 2tt I 

One concludes that the revised version of the thermofield dynamics formulation for fermions introduced in Ref. 
8 plays an important role in order to establish the two-dimensional Fermion-Boson mapping within the thermofield 
dynamics approach. As a byproduct, according with Eq. (|A.15|) , this formulation enable to regard the tilde fields 
of the thermofield dynamics as living on the lower branch of the integration contour in the complex time plane, 
displaced from the real time axis by — i~, in accordance with the work of ref. [7]. 



20 



References 

[1] E. Abdalla, M.C.B. Abdalla and K.D. Rothe, " Non-Perturbative Methods in two dimensional Quantum Field 
Theory", World Scientific, Singapore, 1991, and idem. 2nd edition, 2001. 

[2] L. R. P. G. Amaral, L. V. Belvedere and K. D. Rothe, Annals of Phys. 320 (2005) 399. 

[3] L. Lcplae, F. Mancini and H. Umezawa, Phys. Rep. 10 C (1974), 151; Y. Takahashi and H. Umezawa, Colletive 
Phenomena 2 (1975, 55; H. Matsumoto, Fortschr. Physik 25 (1977), 1. 

[4] H. Umezawa, H. Matsumoto and M. Tachiki, "Thermo Field Dynamics and Condensed States" , Nort-Holland, 
Amsterdam 1982. 

[5] A. Das, "Finite Temperature Field Theory", World Scientific 1997. 
[6] R. Haag, N. M. Hugenholtz and Winnink, Comm. Math. Phys. 5 (1967) 215. 
[7] H. Matsumoto, Y. Nakano and H. Umezawa, J. Math. Phys. 25 (1984), 3076. 
[8] I. Ojima, Ann. phys. 137 (1981) 1; 

[9] H. Matsumoto, I. Ojima and H. Umezawa, Ann. Phys. 152 (1984) 348. 
[10] D. Deleepine, R. Gonzalez and J. Weycrs, Phys. Lett. 419 (1998) 296. 

[11] A. Gomez Nicola and D.A. Steer, Nucl. Phys. B549 (1998) 409; A. Gomez Nicola, R.J. Rivers and DA. Steer, 
Nucl. Phys. B570 (2000) 475. 

[12] A. Liguori, M. Mintchev and L. Pilo, Nucl. Phys. B 569 (2000) 577. 

[13] S. Mandclstam, Phys. Rev. D 11 (1975) 480. 

[14] S. Coleman, Phys Rev. D 11 (1975) 2088. 

[15] K. D. Rothe and J. A. Swieca, Phys. Rev. D 15 (1977) 1675. 

[16] O. Klein, J. Phys. Radium 9 (1938) 1; G. Liiders, Z. Naturforsch. A 13 (1958) 254; H. Araki, J. Math. Phys. 
2 (1961) 267; 



21 



